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Abstract 

The bacterial flagellar motor (BFM) is responsible for driving bacterial locomotion and chemotaxis, fundamen¬ 
tal processes in pathogenesis and biofilm formation. In the BFM, torque is generated at the interface between 
transmembrane proteins (stators) and a rotor. It is well-established that the passage of ions down a transmem¬ 
brane gradient through the stator complex provides the energy needed for torque generation. However, the 
physics involved in this energy conversion remain poorly understood. Here we propose a mechanically specific 
model for torque generation in the BFM. in particular, we identify two fundamental forces involved in torque gen¬ 
eration: electrostatic and steric. We propose that electrostatic forces serve to position the stator, while steric 
forces comprise the actual ‘power stroke’. Specifically, we predict that ion-induced conformational changes about 
a proline 'hinge' residue in a a-helix of the stator are directly responsible for generating the power stroke. Our 
model predictions fit well with recent experiments on a single-stator motor. Furthermore, we propose several 
experiments to elucidate the torque-speed relationship in motors where the number of stators may not be con¬ 
stant. The proposed model provides a mechanical explanation for several fundamental features of the flagellar 
motor, including: torque-speed and speed-ion motive force relationships, backstepping, variation in step sizes, 
and the puzzle of swarming experiments. 

[][]The bacterial flagellar motor (BFM) Is one of only two known protein motors that utilizes the potential energy 
stored In the transmembrane Ion gradient (the ion motive force, or IMF) Instead of ATP, the near-universal cellular 
energy currency. The other such motor Is the one responsible for the synthesis of this energy source—the Fq 
motor of ATP synthase. Understanding how these Ion-driven machines generate useful mechanical work is a 
fundamental issue in cellular biology. Accordingly, much work has been devoted to understanding the torque- 
generation mechanisms of these molecular machines. 

One of the principle diagnostics for a rotary protein motor is the relationship between motor torque and rotational 
speed. Theoretical models attempt to reproduce these empirically measured relationships. The torque-speed 
curve of the BFM was famously shown to display two distinct regimes: a constant-torque plateau at low speeds 
that sharply transitions into a near-linear decrease in torque at high speeds [TJ. However, recent experiments 
from the Berg laboratory show that the number of torque-generating units (or stators) is likely not constant across 
this curve [T?]. That is, stators were shown to come on and offline as the applied load was varied. This is 
akin to a car in which the number of active cylinders changes as the car goes up and downhill. This finding 
effectively invalidates all current theoretical models of torque-generation in the BFM, as they have been based on 
this experimentally measured curve and assumed that the number of working stators is constant (35ll3^[T8lfT7l . 

Recent torque-speed measurements have been performed on single-stator motors [TSl, providing insight into the 
intrinsic torque-generation mechanism of the BFM and making this theoretical reexamination especially timely. 

* These authors contributed equally, 
t Corresponding author: goster@berkeley.edu 


1 



Flagellum 



Figure 1: Schematic showing the basic parts of the BFM. A bacterium has on average four flagellae, each 
attached to the basal body of a motor via a flexible hook. The M-, S-, and C-rings of the basal body are together 
called the ‘rotor’. FUG proteins (26 copies of which are assumed here) are placed around the periphery of the 
C-ring. These interact with the MotA loops of the ‘stator’ to generate torque and rotate the flagella. Stators 
are composed of MotA and MotB subunits, the latter of which attaches the stators to the peptidoglycan layer, 
allowing for torque generation via the MotA-FIIG interaction. A motor can have between 1 and 11 engaged 
stators, depending on the load fRIISOl . 


Furthermore, current models describe torque generation phenomenologically as a free energy surface without 
committing to a specific physical mechanism. In the following we combine the currently available structural infor¬ 
mation with biophysical and biochemical studies on the dynamical behavior of the motor. We use this compiled 
information to propose a mechanically specific and experimentally testable model of torque generation in the 
BFM. The BFM consists of a series of concentric rings embedded in the cell envelope connected to an extracel¬ 
lular helical propeller by a flexible hook (Figure 1). The cytoplasmic C-ring acts as the rotor and the membrane- 
embedded Mot complex(es) act as the stator(s). A motor can have between 1 and 11 such stator units. Each 
stator unit is composed of 4 MotA and 2 MotB helix bundles OUT]. A MotA bundle consists of four membrane- 
embedded a-helices linked by two large cytoplasmic loops. The interaction between these cytoplasmic loops 
with FUG proteins located on the periphery of the C-ring is implicated in torque-generation. We note that while 
there is some controversy in the exact number of FlIGs, this detail does not affect the main points of our model. 
For ease of exposition, in the following we assume that there are 26 FUG ‘spokes’ on the rotor. 

The feat of coupling of an ion gradient to the generation of mechanical work is attributed to the peptidoglycan- 
bound MotB complexes. These complexes each contain an ion-conducting channel with a negatively-charged 
aspartate residue (Asp32), which functions in cation binding. This residue is one of the only strongly conserved 
residues across bacterial species. We propose that the interaction between the Asp32 and a cation passing 
across the inner bacterial membrane (between the periplasm and the cytoplasm) induces conformational changes 
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in the stator complex, resulting in the torque-generating ‘power stroke’ (TTl. While a crystal structure of the stator 
complex is needed for a complete understanding of the power stroke, the available structural knowledge combined 
with information about the motor’s dynamical performance is sufficient to propose a plausible and experimentally 
testable model. Recent structural information suggests that a proline residue in MotA is essential for torque 
generation I4]|TT1[T3]. Using this information, we present a mechanical model for torque generation involving 
proline-induced conformational changes in MotA cytoplasmic loops mm- 

To the best of our knowledge, our model is the first to incorporate known structural information about the BFM 
stator and rotor complexes into a quantitative physical mechanism for the generation of the power stroke. As part 
of this study, we aim to address the following fundamental questions: 

• What are the kinematics and mechanics of the BFM power stroke? 

• What role do charged residues on the stator and rotor play in torque generation, and how does this role 
explain mutational experiments which show only a partial reduction in motor efficiency? 

• What is the shape of the torque-speed curve for motors with a single stator? 

• Why does the motor exhibit backsteps even in the absence of an external ‘reversal’ signal (usually the small 
protein CheY-P)? 

Here we primarily discuss the proton (H+) powered motor of Escherichia coii. However, our model is sufficiently 
general so as to apply to the sodium (Na+) powered motors found in alkalophiles and marine Vibrio species. 


The Mechanochemical Model 


An electro-steric power stroke 

Due to the small magnitude of the forces involved relative to thermal fluctuations, it has long been assumed that 
nearly any form of interaction would be sufficient to explain the rotation of the BFM Q]. For this reason, previous 
models have avoided characterizing the exact nature of these forces, instead treating the interaction between the 
stators and the rotor phenomenologically as a free energy surface and the stator as an ad hoc stochastic stepper 
(33][20]Q9]. 


However, knowledge gained from recent structural (HI^ITg] and biophysical (15] studies have led us to conclude 
that the power stroke of the BFM is eiectrosteric —that is, it is driven by both electrostatic and steric forces. In 
the following, we propose a mechanochemical model consisting of two phases, (i) During ‘electrostatic steering’, 
electrostatic forces position the stator, (ii) Once positioned, the stator delivers a steric push (i.e., a contact force) 
to a FliG protein located along the periphery of the rotor. A more detailed description of the nature of contact 
forces is found in the supplementary text and in reference FTOl . 

In the following we lay out the assumptions involved in the construction of our model, followed by a detailed 
description of the mechanism. Details of the mathematical formulation are provided in the Materials and Methods. 


Electrostatic forces steer the stator into place 


The first step in our model is the steering and positioning of the stator by electrostatic forces. This hypothesis 
originates from the results of the swarming experiments performed by Zhou et ai. |36]. These studies were aimed 
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Figure 2: The predicted energy landscape during electrostatic steering. (A) Schematic of rotor and stator con¬ 
figurations. (j)s and Or are the angular coordinates of the stator and the rotor with respect to the horizontal; 
Ur is the positive angle of the individual FliGs with respect to the radius. Blue arrows denote the direction of 
the dipole (TS] |6]. (B) Predicted surface and contour plots of the electrostatic energy vs. the stator and rotor 
angles. The predicted surface shows the existence of a wide and gently sloping energy well. Note that 4>s and 
Or are periodic variables with periods 7r/2 and 7r/13, respectively; the above plots show one period of each. Our 
calculations consider a single stator centered at (21,-2,1) with the rotor centered at the origin (all distances in 
nm). Computations using this dipole approximation suggest a well of depth --^1 keT for this configuration (see the 
supplementary text for details). 


at elucidating the structure of the MotA loops. They found that mutations of certain charged residues on the 
cytoplasmic portions of the loops hindered—but did not eliminate—motor function. Notably, the deleterious effect 
of mutations on the stator were often able to be countered by corresponding mutations on the FliGs. Certain 
mutations were also found to have very small effects, or even slight improvements, in bacterial motility. 

These results correspond to the idea that mutations of charged residues may result in imperfect steering and 
consequently in a less efficient, but still functioning, power stroke. Likewise, certain mutations may position the 
cytoplasmic loops closer to the adjacent FUG, resulting in a larger power stroke and thus improved motility. 

Note that attractive electrostatic forces strong enough to comprise the entire power stroke would require a non- 
negligible energy to separate the stator and the rotor at the end of the power stroke. This penalty for ‘letting go’ 
would obviate the rotor torque, resulting in a motor with a far lower Stokes efficiency |32] than has been calculated 
for the BFM (~95%) Q]. In contrast, the mechanism we propose here efficiently generates mechanical work from 
the ion-motive force. 

Because detailed structural information of the stator is not yet available, we performed a simple example calcu¬ 
lation to demonstrate how electrostatic interactions can serve to position the stator in place. Details on these 
computations can be found in the supplementary text. For computational convenience, we approximate the im¬ 
portant charge residues on the FUG proteins |6] and stator loops (37] implicated in torque generation as dipoles. 
The assumption that FUG proteins can be modeled as dipoles is based on previous studies (iSl [6]. Modeling 
the electrostatic forces between the stator and rotor by point charge interactions produces results comparable to 
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Figure 3: Ion binding onto Asp32 induces a ‘kink and swivel’ conformational change |[7]. (A) Binding of a proton 
to Asp32 of MotB drives a rapid local reorganization of hydrogen bonds (including those of water). In particular, 
we focus on the creation of a hydrogen bond between the side chain of MotB’s Asp32 and the carbonyl group of 
MotA’s residue 169. Ion binding thus creates a local elastic strain in the MotA helix. The release of this strain 
leads to the proposed conformational change in MotA about the Pro173 residue. Adapted from Kim et al. tTTl . 
(B) Upon ion binding, MotA undergoes a rapid conformational change consisting of three motions: (i) a bending 
about Pro173 </>, (ii) a downward motion, z(^), and (iii) a rotation about its central axis. Inspired by the work 
of Cordes et al. [7], we propose that this ‘kink and swivel’ motion generates the power stroke. Importantly, we 
note that this figure is a two-dimensional depiction of a three-dimensional process, with the motion of the loop 
extending out of the plane of the page. 


those obtained from a dipole approximation. 

Our calculations predict that the energy well produced by the electrostatic interactions will be shallow and wide 
(Figure]^. This results in a weak electrostatic force that is sufficient to position the MotA loop without significantly 
wasting energy to free the stator at the end of the power stroke. Furthermore, the width of the well leads to 
somewhat imprecise positioning. The distribution of observed rotor step-sizes has been shown experimentally 
to be centered around || radians (~13.8°), the average spacing between consecutive FliGs (29] [22]- While this 
result is hardly unexpected, the wide spread of this distribution—in particular, the tendency towards smaller step 
sizes—has been somewhat puzzling. 

Because a wide energy well may result in stators being positioned at ‘nonoptimal’ locations, electrostatic posi¬ 
tioning may contribute to this variance. Since we propose that the stator’s power stroke is imparted via a contact 
force on the rotor, imperfect electrostatic positioning will result in the stator being in contact for only a portion 
of its trajectory. This results in the stator delivering a stroke that is smaller than average. Of course, imperfect 
steering is not likely to be the only factor determining the variance in the observed step-size distribution: the 
uneven spacing of FliG’s along the periphery of the rotor (^, as well as experimental errors, are also likely to 
contribute. 


Motion about a proline hinge provides a steric push 

As proposed previously (H], we hold that the steric portion of the power stroke is the result of a conformational 
change in the cytoplasmic MotA loop. In our model, this motion consists of hinged movements of the MotA 
helices that results in a movement termed ‘kink and swivel’, as shown in Figure]^ (7). The steric mechanism 
proposed below remains valid regardless of which residue, or group of residues, on the MotA/MotB helices acts 
as the inducer. However, we have chosen to focus on MotA’s Pro173 residue because: (i) along with Asp32 on 
MotB, this amino acid is strongly conserved across bacterial species (4], and (ii) previous molecular dynamics 
simulations have found that proline residues induce hinges in transmembrane helices (7], resulting in a movement 
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analogous to the one proposed in the model. The specific mechanism we propose is as follows. 

When a cation binds to the negatively charged Asp32 residue on MotB, the hydrogen bonds (including those of 
water) in the vicinity of Asp32 and Prof 73 on the A3 helix of MotA collectively rearrange. This rearrangement 
induces an elastic strain in the MotA-MotB complex centered around the proline residue in the A3 loop of MotA. 
Figure 1^ shows a candidate scenario, where the carbonyl group of residue 169 on MotA forms a hydrogen 
bond with Asp32 on MotB after proton binding, as proposed in (TTl. This elastic strain induces the kink and 
swivel movement around the proline residue and drives the proposed motion of the lower part of the A3 helix, 
constituting the power stroke (see Figure]^). The binding of the ion and the rearrangement of the hydrogen 
bonds (10^^^ to 10~® s) are near-instantaneous processes as compared to the much slower motion of the kink 
and swivel conformational change (10“"^ s). Thus the chemical steps can be treated as transitions between states 
in a Markov chain. 

The above proposal is supported by a few simple calculations. The maximum torque of the BFM in E. coli is 
~2000 pN-nm jS]. Given that up to 11 torque-generating units may be acting, this corresponds to a maximum 
motor torque of ~200 pN-nm per stator IgQ. As the radius of the motor is --^20 nm, the force generated by a single 
stator during a power stroke is --^10 pN. Direct observation of stepping behavior have shown that the motor takes 
26 elementary steps per revolution, corresponding to a displacement of ~5 nm per step. As explained below, 
our model supposes that each elementary step is actually composed of two half-steps, each imparted by the 
power stroke of a MotA helix. This results in a displacement of ~2.5 nm per power stroke. Molecular dynamics 
studies show the angles subtended by proline hinge motifs from various transmembrane helices to be between 
18-25° (7). From this, we can estimate the length of the cytoplasmic loop measured from the proline hinge to its 
tip to be ~7 nm, a reasonable estimate as the majority of the stator residues have been shown to extend into the 
cytoplasm |[36]. Such a lever arm would result in ~25 pN-nm (--^6-8 keT) of work per half-step, corresponding 
to the rearrangement of 1-2 hydrogen bonds. This energy barrier is sufficient to ensure an efficient directional 
process, as suggested in (9]. 


An in-phase two-cylinder engine 

There are four MotA subunits in each stator complex; see Figure for a schematic of the stator structure. Our 
model supposes that two of these subunits are inactive during torque-generation while the motor is moving pre¬ 
dominantly in a single direction. We base this presumption on the idea that switches between counterclockwise 
(CCW) and clockwise (CW) rotation result from changes in FliG orientation [13l. Given this, we propose that two 
MotA loops are responsible for the power stroke in one direction, while the other two interact with the alternately 
oriented FliG to drive rotation in the other direction. We suppose that loops 1 and 3 are responsible for CCW 
motion and loops 2 and 4 for CW motion, but note that this designation is arbitrary. This mechanism predicts that 
the intrinsic mechanics for power strokes in both directions are equivalent; this has been observed experimentally 


We propose that an elementary step is composed of a pair of power strokes, analogous to the mechanism 
of a two-cylinder engine. It may be possible to observe these substeps in experiments on motors driven at 
extremely low speeds. This can be done using chimeric sodium-driven flagellar motors. As extremes in sodium 
concentration are tolerated far more easily than extremes in pH, these chimeric motors can be driven at very low 
sodium motive forces (SMFs). Thus far, speeds as low as 10 Hz have been obtained l27l . 

A two-ion mechanism can either be ‘in-phase’, in which the energetic profiles of the two stator loops are identical, 
or ‘out-of-phase’, in which their dynamics are offset by a half-cycle. In an experiment using a slowly-driven 
chimeric motor, measuring the rate-limiting step between mechanical substeps can differentiate between these 
two scenarios. For example, if slower ion-binding (e.g., by lowering IMF) extends the dwell time between half¬ 
steps, the out-of-phase engine model is supported. 
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(B) Stator (top view) 



Figure 4: Stator structure and coordinated motion between stator subunits. (A) Proposed arrangement of sta¬ 
tor components as viewed from the periplasm. A stator has four MotA helix bundles, each consisting of four 
a-helices. The four MotA subunits surround a pair of MotB helices. The ion channels associated with the MotB’s 
(shown in green) contain the Asp32 residue essential for proton binding. The stator is attached to the peptido- 
glycan via a linker region on MotB. The power stroke is delivered to the rotor FliGs by the cytoplasmic loops 
between helices A2 and A3 in each MotA bundle (shown as solid bars). Loops 1 and 3 (highlighted in blue and 
red, respectively) are associated with the ion channels. Adapted from Braun et al. and Kim et al. (SlITT]- (B) Due 
to the helical structure of the MotA loops, we can make an analogy between their motion to that of a bundle of 
four gears. Our model proposes that Loops 1 and 3 (shown in blue and red, respectively) drive CCW rotation via 
contact with FUG, while Loops 2 and 4 drive CW rotation. 


The mechanics of these two scenarios are equivalent within the framework of our model. For this reason, we 
discuss only one of these mechanisms in detail: the one in which the two stator loops act in-phase with each other 
(as shown in Figurej^). We choose this alternative because the passage of a single proton across a membrane 
is unlikely to be sufficient to generate an efficient power stroke. A proton passage under standard conditions 
generates ~6 keT, which is slightly less than the calculated length of ‘time’s arrow’ (the energy barrier required 
for a process to be reliably directional in the presence of thermal noise) [9l]. 


A full revolution requires the passing of at least 52 protons 

Our model for torque generation assumes that the rotation of the BFM is ‘tightly coupled’ to the transmembrane 
ion gradient. This means that each elementary power stroke is tied directly to the passage of protons across 
the membrane. Given our prior assumption of 26 elementary steps per revolution, our model thus requires 52 
protons for a full revolution. Previously, a lower bound for the number of ions per full revolution was determined 
by calculating the work done as (r) x 27 r and equating it to the free energy nx IMF, where n is the number of ions 
per revolution and IMF is the ion motive force, as before [15]. The above calculation resulted in an estimate of 
n = 37, lower than the 52 ions per revolution supposed by our model. 

This discrepancy can be explained as follows. While the above is indeed a lower bound, a tighter bojjnd can be 
computed. The calculation of work as stated above suggests that the power output per revolution is P = {t){uj). 
However, power is formally calculated as P = {t ■ w), which differs from P by a covariance term, cov (r, w). This 
follows from the fact that, for any two stochastic processes X and Y, (XY) = (X) (Y) + cov{X, Y). 

Note that the number of protons per revolution assumed by our model is also a lower bound; that is, we have 
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assumed that 52 working ions are required per revolution. Many factors can result in the passing of more ions 
than predicted, including leakiness of the ion channels, irregular arrangement of FliG’s around the rotor, or im¬ 
perfect placement of stators by the electrostatic steering forces. This can be quite easily extended within our 
mathematical framework by replacing the step function associated with ion binding with a sigmoidal function. 


The mechanical escapement 

Figure!^ depicts the mechanics associated with the power-stroke. We choose the angle subtended by a stator 
loop (j)g (where i corresponds to the loop number) as the order parameter. That is, we consider the energy 
landscape along the arc length of the mechanical trajectory of the stator loop. A stator loop has two stable 
configurations: straight {(pg = 0°) and bent {(j)g 20°). Both of these configurations correspond to energy 

minima in different chemical environments: when the negative Asp32 is not neutralized by a proton, the loops 
prefer to maintain a straight posture {(j>g = 0°). The presence of bound protons induce a free energy change 
sufficient such that a thermal fluctuation can induce the conformational change to the bent state {4>s ^ 20°). 
During a power stroke, the entire stator complex undergoes a collective gear-like motion as shown in Figure]^. 
The conformational change due to the ‘hopping on’ of the ion produces the first half of the power stroke: here, 
loop 1 pushes the FUG, while loop 3 is put in place to carry out the second half of the power stroke during the 
reset (Figurej^). This ‘reset’ corresponds to the ‘hopping off’ of the proton, resulting once again in the the stator 
loops surmounting the energy barrier between configurations and reverting to the straight position (</)! 3 = 0°). 
Note that the numbering of the loops is arbitrary; the mechanism proposed here is equivalent to one in which 
loop 3 performs the first half of the power stroke and loop 1 the second. 

In summary, a torque-generation cycle by a single stator of the bacterial flagellar motor proceeds as follows: 


1. Electrostatic interactions between charged residues on MotA and FUG steer a stator tip close to a rotor 
FUG. 

2. In the presence of a membrane potential, the two MotB aqueous ion channels open and two protons bind 
to the negatively charged Asp32 residues on the MotBs. This triggers a reorganization of the hydrogen 
bonds in the vicinity of the Pro173 on MotA (see Figure]^). 

3. The hydrogen bond rearrangements induce elastic strain in the straight MotA loops. This strain drives a 
‘kink and swivel’ motion of the MotA loop, increasing the bend angle (from = 0° to 20°, as shown in 
Figure|^). 

4. One MotA loop (Loop 1, shown in blue in Figurej^) applies a steric push to the nearest FUG, resulting in 
one-half power stroke. 

5. At the same time, the movement of the stator ion binding pocket moves downward so that the pocket is 
exposed to the cytoplasm. The ion-channel is now closed to the periplasm. The protons ‘hop off’ MotB into 
the cytoplasm, now inverting the strain in the bent MotA loops . 

6. The inverse strain drives the movement of the loops in the reverse direction, straightening the bent MotA’s 
(i.e., from (j)g ~ 20° to 0°). 

7. The other MotA loop (Loop 3, shown in red in Figure [^) now applies a steric push to the same FUG, 
completing the second half of the power stroke. 


Consequently, to the rotor, the stator appears to be an ‘inch-worm’ stepper with FliGs as the ‘stepping-stones’. 
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Figure 5: Dynamics of the rotor-stator interaction. (A) Mechanics of the power stroke. Top panel: After 
the initial electrostatic steering, two protons bind to the charged Asp32 residues on the MotB’s. The consequent 
rearrangement of hydrogen bonds induces an elastic strain in the straight MotA loops. Release of this strain 
results in synchronous ‘kink and swivel’ motions about the proline ‘hinge’ in both MotA’s. As a result, a steric 
push is imposed on FUG, and the first half of the power stroke performed by Loop 1. Importantly, this motion also 
has a vertical component—the loops lower themselves out of the membrane. Bottom panel: The lowering of the 
MotA loops exposes the protons in MotB to the cytoplasm, where they are released. This results in a ‘reset’ of 
the MotA loops, during which the Loop 3 carries out the second half of the power stroke. We note that this image 
depicts a two-dimensional projection of a three-dimensional motion: the motion of the stators is not constrained 
to the plane of the page. (B) Energetics of the power stroke. Because the two loops move in-phase with 
each other in our model, their energetic pictures are identical. We describe the free energy landscapes using 
double-well Landau potentials. These landscapes are shown in blue for Loop 1 and red for Loop 3 with respect 
to the angles of the stator 4>, and rotor 9. The initial entrance of the proton into the ion channel (kon) places the 
system within ksT of the energy barrier. Thermal motions then result in the first half of the power stroke {Top 
and middle panels). The exit of the protons into the cytoplasm (fcotf) results in the ‘reset’, and the second half of 
the power stroke {Middle and bottom panels). 


Results and Predictions 


Using the mathematical model described in the Materials and Methods, we performed both analytic calculations 
and numeric stochastic simulations. Statistics from simulated trajectories—an example of which is shown in 
Figure —were used to calculate various experimental quantities including average motor torque and angular 
speed. 

In the sample trajectory for the rotor motion, the duration of a power stroke (T^) and the waiting time between 
consecutive power strokes (r^,) are highlighted in orange and purple, respectively. The highlighted power stroke 
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Figure 6: Summary of recent experiments and comparisons with simulations of our model. Results in (A,B,C) 
are derived from numerical simulations. Results in (D) are analytic curves obtained from a deterministic approx¬ 
imation (see the supplementary text). In all plots, model calculations are shown by solid lines, and experimental 
data is shown as open colored circles. (A) Sample trajectory generated by the model. Moving (Tm) and waiting 
(Tn/) times are shown with orange and purple backgrounds, respectively. Occasionally, reversals (shown with 
green background) appear when MotA Loops 2 and 4 are engaged due to conformational changes in FUG. (B) 
Single-stator torque-speed curves measured in a chimeric sodium motor for various sodium concentrations at 
pH 7.0. Curves show a concave-down shape, with the length of their ’plateau’ being SMF-dependent (data from 
Lo et at. [15]). (C) Motor speed vs. SMF in a chimeric sodium motor shows a nearly linear relationship across 
various loads (data from Lo et at. (T5]). (D) The maximum speed at low and high load may show different rela¬ 
tionships with the number of stators (data from Sowa and Berry [28]). We note that this result is dependent on 
the assumption that all stators work in synchrony, and that it has come into question given recent results from the 
Berg laboratory lfT4l . 
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shows two half-steps, corresponding to the two sequential steric pushes by the two MotA loops involved. As in 
experimental trajectories, occasional reverse steps are also observed in our simulations, one of which is shown in 
the sample path in Figurej^. An explanation for back-steps that is compatible with our model is provided below. 

The results shown in Figure[^,B,C were obtained via simulation. Analytic calculations on an approximate deter¬ 
ministic model (explicitly provided in the supplementary text) were also performed. The results shown in Figure 
were computed analytically. Finally, in addition to fitting current experimental data, we propose further experi¬ 
ments to test, and hopefully validate, several of our model predictions (e.g.. Figure]^. 


Single-stator motors exhibit concave-down torque-speed curves 

Until recently, BFM experiments have been performed on motors with multiple stators, with no direct accounting 
for the number of engaged stators at a given load. Therefore, the existence of the torque-speed ‘plateau’ and 
‘knee’ have been assumed to be innate characteristics of the rotor-stator interaction, largely because there was 
no evidence to the contrary. However, Lo et at. [TS) performed experiments using a chimeric single-stator motor 
showing smoother torque-speed curves without a dramatic plateau as observed for wild type motors. While these 
curves are still concave-down in shape, the extent of the plateau regions are quite variable, and are dependent 
on the IMF. 

Our simulations show torque-speed relationships consistent with these single-stator experiments (Figure |^). 
The behavior of the torque-speed curves results from a competition between the time taken for a mechanical 
half-step (Tm) and the waiting time between ion-binding events {Tw)- For example, our simulations show that 
the average time in moving a half-step (Tm) is ~2 to 20 ms at high loads and --^0.01 ms at low loads. The 
average waiting time under standard conditions {Tw) is --^0.2 ms |20]. Therefore, at low loads, the motor is in 
a kinetically limited regime, where the waiting time between steps is generally higher than the time required to 
complete a step. Conversely, the motor is mechanically limited at high loads when (Tm) > {Tw)^ resulting in the 
observed plateau. Consequently, as shown in Figurej^, this plateau region grows smaller as the IMF decreases 
(i.e., as {Tw) increases). 

This competition is also manifested in the relationship between speed and IMF: speed depends linearly on IMF at 
high loads, but in a slightly nonlinear fashion at low loads (Figure]^). Given that the rotor moves 27r/26 radians 
per step, the speed of the rotor {ujr) can be approximated as: 

_ Stt 1 

^ 26 (Tm) + {Tw) 

At high loads, (Tm) » {Tw)- Because the force applied during a power stroke is inversely proportional to the 
ion-motive force, {ujr) cx IMF at high loads. In contrast, the waiting time eclipses the time for a mechanical step 
at low loads, and therefore ujr cx l/{Tw) oc exp(IMF/fcBT). Further details to this end are provided in the 
supplementary text. 


Backstepping in the absence of CheY-P is due to thermal flipping of FliG 

The BFM plays a central role in bacterial chemotaxis: the direction of rotation of the motor determines whether a 
bacterium will move in a straight line (CCW) or ‘tumble’ (CW) to move in a random new direction. This switching 
is typically initiated via a signal transduction pathway, in which a response regulator protein, CheY, is phospho- 
rylated into an activated form, CheY-P, to induce tumbling. For more information on this pathway and bacterial 
chemotaxis, we refer the reader to several excellent reviews 13111241 . 
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However, occasional backstops (e.g., CW motion during primarily CCW rotation) are observed even in the ab¬ 
sence of CheY-P. This has been attributed to microscopic reversibility, of which three possible models are dis¬ 
cussed in the supplementary text. For example, Mora et al. ascribed switching in the BFM to the diffusive motion 
of the rotor through a ‘bumpy’ 26-fold periodic potential [21]. However, recent structural studies have found that 
there exist two main configurations for the FliGs (T^ [22], lending support to the idea that a ‘flipping’ between 
these states is the molecular basis for backstepping. We note that despite a general agreement on the existence 
of two distinct FliG configurations, the exact nature of the conformational change to the CCW direction remains 
controversial. 

In our model, the probability of observing a backward step is equivalent to the probability of finding a FliG oriented 
in the CW state (assuming a primarily CCW-rotating motor). Within the framework of our model, whenever a FliG 
changes its state and is close to a stator, then the stator uses loops 2 and 4 to apply a contact force and pushes 
the FliG in the CW direction. To model the flipping between CW and CCW states for the FliGs, we use a nearest- 
neighbor periodic Ising model with the 26 FliGs arranged on a one-dimensional ring. Such models have been 
used successfully to explain rotational switching (see, e.g., O). 

In our model, when the FliGs are oriented at an angle of roughly 10-20° with respect to the radial direction, as 
shown in Figure]^, the motor moves in the CW direction by virtue of contact forces from loops 2 and 4. Con¬ 
versely, when the FliGs are pointed either orthogonal or at an angle of 180° with respect to the CW orientation, 
the motor steps in the CCW direction using loops 1 and 3. Using an Ising model for the flipping of FliGs, we 
calculate the probability of a backstep to be ~8%. This probability was calculated to be ~7.3% from stepping 
statistics collected by Sowa et al. demonstrating that a backstep might indeed be simply due to fluctuations 
in FliG orientation. Further details on these calculations are provided in supplementary text SI. 


Speed and stator number are nonlinearly related at both high and low loads 

Recent experiments have shown that the bacterial flagellar motor is a dynamic structure, constantly exchanging 
stators in response to changing viscous loads (I4][l2l[^. Recently Tipping et al. ||30] showed that the number 
of torque-generating stator units depends on the viscous load. Dynamic studies have suggested that new stators 
may come online as the load increases P3- This result appears to be responsible for the ‘resurrection’ experi¬ 
ments which show a stepwise increase in speed as stator units are recruited at high loads, but only a single jump 
at low loads (28]. 

Under the assumption that individual stator units act in synchrony with each other, the mathematical model pro¬ 
vided in the Materials and Methods can be extended to deal with the dynamics of motors with multiple stators. 
The details of this extension and generalized equations for motors with N stators are provided in the supple¬ 
mentary text. The speed shows a slightly nonlinear dependence on the number of stators at very high loads 
(Figure |^). A linear relationship between speed and stator number would imply that the applied force of the 
stator is independent of the load; our results suggest that this nonlinearity may arise as a natural consequence 
of the steric force. 

Furthermore, our computations match results from the resurrection experiments very well, even when a constant 
number of stators is assumed (Figure]^). Given the assumptions of the extended model, we predict that the 
motor speed at low loads will be independent of stator number if (and likely, only if) individual stators step in 
synchrony. Further experiments are needed to determine if this is indeed the case. Because of the recent 
experiments performed on single-stator motors, we have chosen to primarily focus on the innate mechanism of 
torque-generation in the BFM. For this reason, a more thorough investigation of the interaction between stators, 
while worthwhile, is left for a future study. 
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Rotation rate, CO [Hz] 


Figure 7: Prediction of how the number of activated stators at a given load and torque-speed relationship in a 
constant-stator motor can reproduce observed torque-speed curves. Red open circles show data from Yuan and 
Berg for motors movithe supplementary textng in the (A) counter-clockwise (CCW) and (B) clockwise (CW) 
directions. Data was collected from motors in which the number of engaged stators was not constant across 
loads. Analytic curves from our model for motors with N = 1,..., 7 stators are shown in blue. The asymmetry 
between the observed CCW and CW torque-speed relationship is likely due to differences in stator recruitment 
PQ]. 


The observed constant-torque plateau is due to stator recruitment 


Dynamic experiments by Lele et al. showed that the number of active torque-generating units in the BFM is not 
constant over all applied loads (14], thus bringing into question current interpretations of the well-known torque- 
speed curve (^. Our model simulations are consistent with the single-stator curves measured by Lo et al., 
lacking the large constant-torque plateau seen in previous experiments where the number of stators was not 
directly accounted for. 

Given this, we suggest that the observed two-regime curve is a consequence of variable stator number rather 
than the rotor-stator interaction. Torque-speed curves predicted by our extended model for varying number of 
stators are shown as solid blue lines for iV = 1,2,..., 7 stators, overlaid with experimental data (red circles) 
for CCW (Figure [^) and CW (Figure |^) motors. The relative motor torque is shown as in the experiment: 
normalization is performed relative to the torque generated by a motor with 6 stators moving at ~70 Hz (the 
lowest speed measured in for the CCW motor). 

In particular, as shown in Figure]^, we predict that the ‘constant torque plateau’ seen in CCW motors is due 
to the motor maintaining a constant number of stators, and the sharp transition to decreasing torque is due 
to stators becoming inactive at a critical load or speed. Furthermore, as shown in red in Figure |^, torque- 
speed measurements on CW motors yield a near-linear torque-speed relationship (^. In accordance with our 
prediction, this asymmetry was reflected in the experiments of Tipping et al.: CCW motors cease to recruit 
additional stators at loads far lower than CW motors |^. This suggests that CW motors ‘jump’ between curves 
as they recruit additional stators at low loads, but stay on a single curve (i.e., maintain a constant number of 
stators) once the load is sufficiently high, while CW motors continually recruit stators at loads far lower than CCW 
motors. 

To this end, our results support the idea that the characteristic ‘knee’ commonly associated with the dynamics of 
the flagellar motor may simply be an artifact of experiments performed on motors with a variable number of stator 
units d?], rather than an intrinsic feature of the rotor-stator interaction. This further asserts that a model for the 
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loading and unloading of torque-generating units is needed to completely interpret the torque-speed curves for 
experiments where the number of active stators was not monitored. 


Discussion 


The ability to convert a transmembrane ion gradient into rotary torque is rare, observed in only two known protein 
motors: the Fq motor of ATP synthase and the BFM. The mechanism behind the torque generation in the latter 
has been a longstanding mystery, driven by the fundamental role of this machine in bacterial locomotion and 
chemotaxis. 

Recent experiments have shown that the number of torque-generating units in the flagellar motor is likely dynamic 
and not constant across previously measured torque-speed curves (14]. This result brings into question previous 
models of torque-generation that have relied on these curves. Here we combined recent biophysical [15] and 
structural [HlQl] studies to reexamine these models. 

Crucially, we have utilized all known structural information of the BFM, as well as the experimental measurements 
on single-stator motors by Lo et at. [15], to present the first mechanically specific model of torque generation. 
Using this information, we are able to present an explicit model of the dynamics of the stator during a torque- 
generation cycle. Our model implicates a steric interaction between the cytoplasmic MotA stator loops and the 
FliG proteins of the rotor. We have predicted that this interaction is driven by conformational changes in the 
MotA loops due to the binding of a cation to an essential aspartate residue on MotB. Results from our model 
simulations reproduce recently measured torque-speed and speed-IMF curves from single stator motors. 

The model we have proposed is akin to a two-cylinder engine, where two of the four MotA loops act when the 
motor is moving in the counterclockwise direction and the other two loops act in the clockwise direction. We have 
proposed that the two loops act ‘in-phase’ with each other, moving in synchrony as two protons bind to the MotBs 
and are subsequently released into the cytoplasm. In this manner, the first loop executes its half of the power 
stroke when the protons are bound to the MotBs, and the second does so once they have ‘hopped off’ into the 
cytoplasm. 

Experiments performed at low IMF can be used to differentiate between a stator acting as an engine which is 
‘in-phase’ or ’out-of-phase’. Because ion-binding is rate limiting under these conditions, trajectories would show 
clear half-steps only if the BFM acts as an ‘out-of-phase’ engine. However, given that the mechanics of the power 
stroke for both scenarios are equivalent, the mechanism corresponding to an out-of-phase engine would lead to 
a calculation analogous to the one presented in this work. 

A notable prediction of our model arises from extending our equations for single-stator motors to model the 
dynamics of motors with multiple stator torque-generating units. Our results predict that the observed two-regime 
torque-speed curve is likely the result of a load-dependent recruitment of stators (Figure [^), with the sharp 
transition between regimes corresponding to the point at which stators begin to come off the motor. This lends 
strong support to the theory that the unique shape of the BFM torque-speed relationship may only exist due 
to the inability of most experiments (until now) to explicitly monitor the number of engaged stators at any given 
load |14l- This further indicates that a model of stator recruitment is necessary to fully understand experimental 
results performed on motors where the number of stators was not necessarily constant. Furthermore, it suggests 
that low-load measurements were never successfully performed on motors with multiple stators engaged, thus 
shedding doubt on the result that increasing stator number does not effect motor speed at low loads. 
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Table 1: List of parameters used with units, values, and reference. 


Parameter 

Definition 

Units 

Values 

Ref 

R 

Radius of the rotor 

nm 

20 

niT 

£p 

Length of the proline hinge arm 

nm 

7 

m 

Cs 

Drag coefficient of the stator 

pN-nm-s-rad~^ 

0.0002 

- 

Cr 

Drag coefficient of the rotor 

pN-nm-s-rad~^ 

0.017 

- 

Cl 

Drag coefficient of the load 

pN-nm-s-rad~^ 

0 .0-10 

ED 

(t>s 

Angular position of the stator 

rad 

- 

- 

0R 

Angular position of the rotor 

rad 

- 

- 

Ol 

Angular position of the load 

rad 

- 

- 

K, 

Hook spring constant 

pN-nm-rad^^ 

1000 

El 

N 

Number of stators 

- 

1-11 

m 

T 

Rotor torque from stator 

pN-nm 

- 

- 


Electrostatic potential 

- 

1.5-2 keT 

- 


Methods 


The mechanochemistry of the torque-generation cycle of a flagellar motor with a single stator unit can be modeled 
by the following Langevin equations. The dynamics of the angular positions of the stator loops * G {1 j 3} 
are given as; 


Cs 


dt 


dG dVRs 

Torque from Reaction 
Proiine hinge from rotor 


+V2kBT(sfnit). 
o(l>s " ---" 

V—^ Thermal fluctuations 

Electrostatic 
attraction 


( 1 ) 


Here, as in the following equations, the last term is the stochastic Brownian force, where N{t) is uncorrelated 
white noise. Cs is the effective drag coefficient of the stator. G = G{(j)g,j) denotes the free energy of stator 
loop i, modeled in Figure as a Landau potential. However, because of thermal fluctuations the exact shape 
of the potentials is immaterial. Accordingly, we approximate this potential by a piecewise linear function for ease 
of computation. The parameter j g {0,2} corresponds to the chemical state of the system: j = 2 if two protons 
are bound to the MotB helices and j = 0 if not. The switching between the two chemical states corresponds to a 
‘jump’ between potential curves, as shown in Figurej^. 


As the stator moves between the two configurations, it induces a contact force, and subsequent torque, on the 
rotor. Unlike previous models, we do not assume that this torque is constant across loads, but rather depends on 
the Cl (see the supplementary text for more information). To this end, we do not allow a linear interaction potential 
between the stator and the FliG; this would result in a constant applied force, which is not true for contact forces. 
We model the steric interaction potential Vrs as 


Vrs{4>S^ ^r) 



(ROR—tptp's) 

Xrs 


if 0 < X < Xrs 

otherwise, 


where x = Xrs -F R0r - ipcjig denotes the distance between the position of the stator loop and the nearest 
FliG. From this, the torque imposed on the rotor is calculated as Tcontact = -dVRs/dOR, while the corresponding 
reaction torque on a stator loop is given by Treaction = -dVRs/d(j)s- 

The charges on the FliG and the stator loop exert weak attractive forces on each other. These forces prevent 
the drifting of the rotor with respect to the stator during the chemical transition events. We refer the reader 
to the supplementary text for more on the effects of attractive electrostatic forces on the torque-speed curves. 
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With the contact torque and the weak electrostatic forces, the total instantaneous torque on the rotor is given by 
T = "rcontact - R{d'ip/d9ii)- The average torque on the rotor is a time (or ensemble) average of the instantaneous 
torque. Finally, the rotor and load are connected by a linear spring with elastic constant k; the elastic coupling 
terms in the equations for the rotor and the load thus appear with opposite signs. 

Given this, the rotor dynamics are described by a corresponding Langevin equation: 

+ ( 2 ) 

ot/n oon ■. -^-/ y--/ 

' ^ ^ Connection Thermal fluctuations 

Torque E.S. to load 

from stator attraction 

where Ca is the effective rotor drag coefficient. Finally, the dynamics of the load are then driven by the motion of 
the rotor: 

Cl^ = KiOn - Ol) + y/2kBTaUt).^ (3) 

Spring connection Thermal fluctuations 
to rotor 

As above. Cl is the effective drag coefficient of the load. 

The above model can be collapsed to explicitly include only the dynamics of a stator with a single loop that 
generates torque both during its bending (0^ increasing) and unbending (^>5 decreasing). This description is 
isomorphic to the mechanism described previously (see Figure]^ because the mechanics of the two half power 
strokes are equivalent as described above. The equations corresponding to this reduced model are provided in 
the supplementary text. 


dOa _ 
dt 
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1 Model equations 


In the Materials and Methods section in the main text, we provided Langevin equations describing the dynamics 
of the two stator loops, the rotor, and the load. We explicitly model the dynamics of the stator loops, the rotor, and 
the load. As described previously, the mechanics of the two stator loops are equivalent. Therefore, a reduced 
model describing the dynamics of a stator with a single loop that generates torque both during bending (^s 
increasing) and unbending (t/fs decreasing) is isomorphic to the one described in the main text. 


Stator: 


Rotor: 


Load : 


. d<j)s 


Cr 


R 


d0 

dt 


Cl 


d^L 

dt 


dG{^s,j) , _ ^ 

dc^s 5</>s d(^s ^ 

Torque from Reaction Electrostatic 

Proline hinge rotor attraction 


dVRS 

Torque 
from stator 



Electrostatic 

attraction 


Rj^R — ^ l ) 

Spring connection 
to load 


RiOR-0L) + V2kBT (: Lfnit)._ 

Spring connection Thermal fluctuations 
to rotor 


yj2kBnsUt) 

Thermal fluctuations 


V2kBTCnfnit) 

'-V-' 

Thermal fluctuations 


( 1 ) 

( 2 ) 

(3) 


Here, as before, (s, Cr^ ^^d Cl are the effective drag coefficients of the stator, rotor, and load. The last term in 
each equation is the stochastic Brownian force, where /„(t) is uncorrelated white noise. All other symbols are as 
described in Table 1 in the main text. In writing the equations of motion, we assume that the stator rotates about 
the proline hinge. Similarly, the rotor rotates about the axis normal to the plane of the rotor and passing through 
its center. 


In Equation fi], the internal force driving the stator due to the rearrangement of hydrogen bonds caused by a 
proton binding event is denoted by Fp = Here, G{(j)s,j) denotes the free energy of the stator. Due to the 

fact that thermal fluctuations are of a comparable magnitude to the free energies considered, the exact shape of 
the potentials is relatively unimportant. For ease of computation, we approximate the potential using a piecewise 
linear function. In this setup, the force applied on the rotor by the stator loop is constant and positive during 
each mechanical power stroke. At other times, there is little elastic strain on the MotA loops, and accordingly, the 
applied force is near zero. 


The torque generated by the stator is dependent on the applied load (l as a natural consequence of steric 
forces. A general discussion on contact forces and the explicit formulation of the repulsive interaction potential 
Vrs, are provided in Sectionj^and Section [TT!2l respectively. The contact torque applied to the rotor (in Equation 
l^), and consequent reaction torque applied to the stator (in Equation |[i]), are given by rcontact = ^ad 

_ _ _ OVrs 

'reaction — ■ 

Additionally, the charges on the FUG and the stator loop also exert weak attractive forces on each other. These 
forces prevent the drifting of the rotor with respect to the stator during the chemical transition events. Further 
discussion of these forces are included later in Section]^ In the above equations, we denote this term via the 
term --^R; it detracts from the repulsive torque imposed by the steric force of the stator. 


The rotor and load are connected by a linear spring with constant k; the elastic coupling terms in the equations for 
the rotor and the load thus appear with opposite signs (in Equations ||^ and ||^, respectively). The elastic constant 
in the experiments can vary depending on the length of the hook when attaching the bead. In some cases, the 
hook is very short or is stiffened by an antibody linker. This corresponds to a large spring coefficient [?,0. An 
analysis of this model in the corresponding limit k oo is provided in Section pi~T2|of this supplementary text. 
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1.1 Simplified deterministic model 


By taking an average over many trajectories, it is possible to generate a deterministic analogue of the model 
presented above. Although numerical simulations on the full stochastic model were used for the results in this 
manuscript, the below formulation is convenient primarily for expository purposes. In particular, it admits explicit 
analytic solutions for many experimental situations. Before we provide the numerical implementation of the full 
Langevin equations Q-j^, we use the following model to introduce several important concepts. 


The deterministic equations of motion can be obtained by time averaging the equations as 

C,S ^ tp 

= {t) - K{eR-9L) 


(4) 

(5) 

( 6 ) 


Note that in addition to time averaging, we have neglected the electrostatic term for computational convenience, 
as it tends to be quite small in value. Here, the average torque on the rotor (r) results from averaging the torque 
on the rotor as Tcontact = many trajectories. The return force then can be calculated by (r) /R, which 

is then multiplied by (.p to calculate the return torque. The internal torque of the proline hinge for a given ion 

dG 

motive force (IMF) is — oc IMF. The averaged equations do not contain a noise term because the 

d(j)s 

terms /„(i) are Gaussian with mean zero. In the following, we compute expressions for the average torque and 
speed during a single power stroke of the motor from the above deterministic model. 


Under the assumption that all active stators act in synchrony. Equation Q can be generalized to a motor with N 
stators as follows. An analogue of Equation 0 now corresponds to the motion of the i-th stator (i e 1,2,... N)\ 


Cs^ = Fpip-^^ep, l,2,...,A^. 


We then can sum the equations of all stators, which results in 

dcj)]. 


NCs- 


dt 


= NFpip - 


R 


(7) 


( 8 ) 


Note that we have used the fact that N{Ti) = (r). Similarly, an equation for the rotor in a motor with multiple 
stators can be written as: 


dOR 

dt 


N 

Y^{F)-R{en-0p^) 

2=1 


(t) - K( 6 »i^ - 9l), 


(9) 


which is the same form as Equation ||^. Note that only the terms corresponding to the stators are summed (i.e., 
the connection term between the rotor and the load is not multiplied by N). Equation |6| also remains as for the 
single stator case because the spring connection term is unaffected by the addition of torque-generating units. 


The three equations of motion ID-lD for a single stator contain 4 unknowns {9r, 9l, (j)s, (r)}. This results in an 
indeterminate system, and requires the addition of an equation to generate a unique solution. This additional 
equation can be obtained from a fundamental property of contact forces. Since the stator loop is in contact with 
the rotor during the power stroke, the velocities of the stator loop and the rotor must be equal. This leads to a 
‘contact condition’ for the tangential velocities of the stator loop and the rotor: 


, d<h _ j^d9R 
P dt ~ df 


( 10 ) 
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1.1.1 Analysis 


Multiplying Equation ||^ by ^ and summing with gives: 


Rearranging we get 
Substituting into equation fiT) : 


NRCsd(j)s , ^ d0ii fa a \ 


d(t)s _ R dOft 
dt £p dt 


NB?C,sd0R dOfi I a a \ 

-If- ^ 

dOR (NR^Cs , -1^ ATT? v (a a \ 

+ V = 

CR^iNM + 1 ) = NFpR - Ki0R - 9 l), 


where M = Rewriting gives us 


dO 


R 


NFpR 


dt Cr{NM + 1) CrINM + I) 

Dividing equation by Cl and subtracting leads to 


didR - Or) 


FpR 


dt Cr{NMfI) 

This differential equation has the solution: 

{Sr - dL){t) = 


1 


{Or-Or). 


1 


Cii(iVM+l) Cl 


+ {Or - Or). 


a’ 


where A* := 


Cr(nm+i) 


and a* := 


(cii(WM+l) + Ct)- 


( 11 ) 


( 12 ) 


(13) 


(14) 


(15) 


(16) 


Substituting expressions for and from Equations ||^ and respectively, into the equation for contact 
condition |[T^ leads to: 


dt 


R 


{NFpR - r) = — [t - k{9r - 9r)]. 


'' Cr 


(17) 


Note that, because we are interested in the full time course, we consider r = T{t), rather than the average torque 
(r). Rearranging and solving for r: 


NCrR [Fpfp - CsRr{0r - Or)] 

^ “ -(Rf+NCsI^ ■ 

Plugging in the derived expression for {9r - Or) from Equation p6) : 

CR{tp(C,L+CR) + NR?Cs) 


(18) 


FpifRCRCs 


Cs 


£Kr + iVi?2Cs 


(19) 
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Figure 8: (A) Rotor speed and (B) Torque as a function of the load drag. Cl- Plotted from Equations and 
| |20) , respectively. 


From this value, we can calculate the average torque in a single step (r) by 


1 /■^- 

(r) = ^ / T 

JQ 

FptlNR ( 

{t)dt 

^-1 + exp 

KTm 

( Cl CR+NR^s/el) 

]) 

NR^lCs + kT^{Cl + Cr) KiCL + Cr) + NR^s) 





{eiiCL + Cr) + NR^Cs)" 


( 20 ) 


dO]^ 

where Is the time spent moving during a step. Likewise, we can also calculate the speed of the load from 
Equations and ffS) : 


dt 




1 — exp 


— Kt ^ + 


a ^ Cr+NRHs/^: 


R\s/ei) ) 


FpilNR 


eua + Cr) + NR^Cs 


( 21 ) 


As with torque, we integrate over a time step to find the average speed 




/ d9L\ FpipNR 

(^-1 + exp 


( Cl Cr+NRHs/Ci) 

]) 

Cl {£ICr + NR\s) + {L\Cl + Cr) + NR\s) 


\ dt / 

KTm{£liCL + CR) + NR^Csf 


( 22 ) 
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Using Equations ||^ and we can calculate a family of parametric torque-speed curves (parametrized by the 
load Cl), where each curve corresponds to a motor with a constant number of synchronously stepping stators. 
Computed curves for motors with 1,2,..., 7 stators are shown in Figure 7 in the main text. 


1.1.2 Approximation: modei without spring 


In most experimental setups, the filament is removed and a bead is attached to a shortened hook connection. 
Additionally, the hook is sometimes stiffened with an antibody linker. These setups have a rigid connection 
between the rotor and the load, corresponding to a large spring constant. In this section, we perform a similar 
analysis to that in the previous section for the limit k oo. The calculations performed in this section provide 
analytic formulas for a clear physical understanding of several important properties of the model. 


The rotation rates of the rotor and load become equal after an initial “wind-up” period. That is the rotor and load 

/ dO dOi, \ 

move at the same angular speed (i.e., after the system reaches a steady state. Note however 

that the angular positions 9^ and 
Equations and 


dt dt 7 

L Still maintain a (constant) offset. This can be seen explicitly by subtracting 


d6ii d6L 

dt dt 
djOR - Ol) 
dt 


(r 

Cr 


1 

Cr 

1 

Cr 


1 

Cl 

1 

Cl 


{9r - 9l) 
{9r - 9l)- 


(23) 


To simplify some notation, we define x := {Or - Or), a := k 


j_ 

Cl 


and A ■= 


and rewrite Equation 


dx 

dt 


= A — ax. 


As in the previous section, we can solve for the timecourse x{t) = {9r - 9L){t): 


(24) 


x{t) = {0R-eL){t) = ^^^ ^ 


(25) 


By definition, a > 0, and so it is clear that x{t) will reach a constant value after an initial startup. As k ^ oo, the 
‘wind-up’ time goes to zero and ^ ^ in this limit. 


Summing Equations and ||^ gives us: 


dOR 

dt 


+ Cl 


d^ 

dt 


(t)- 


(26) 


We approximate ^ ^ when k is large, and Equation reduces to: 


{CR + CL)^ = {r). 

We multiply Equations and by £p and R, respectively. After some algebra: 


, dC’s 
^ dt 
dOR 


£2 


NF — — 

NCs \ ^ R 


R 


R{r) 


dt 


(Cl + Cr) 


(27) 


(28) 

(29) 
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Figure 9: Angular difference between rotor and load. (A) The change in {On - Ol) with time decreases steadily 
until a fixed point is reached (marked in red). (B) Accordingly, {On - On) increases and asymptotically reaches a 
steady state. The steady state is denoted by a red line, corresponding to the zero of shown as the red 

dot in (A). 


Given the relationship {T^, we equate the two right hand sides of and j2^ : 


NCs 


NF,-^\ = 


R{t) 


{Cl + Cr) 
NR^{t) 




(30) 


Here, the second line is obtained simply by multiplying through by N and R. Solving for (r): 



NFpR 


v) — 



' {Cl+ChKJ 


(31) 


We can use this expression to attempt some intuition for the result of Sowa et al. regarding torque and stator 
number at high and low loads. Consider the following two limits of Equation (i) high load, when Cl > Ci?> 
and (ii) low load, where Cl < Cr- 


In the first case, we have: 

('r)high = 


« NFpR. (32) 

The third line follows from the fact that Cs < Cr ^ Cl- For very high loads, the observed torque is 180 pN- 
nm, and therefore the force Fp « 9.5 pN. Equation suggests that the torque increases linearly with stator 
number under extremely large loads as observed in the experiments. However, for a given Cl, it can be seen 
that nonlinearities can arise in the torque versus number of stators even in the high load limit. This property is 
primarily due to the nature of contact forces, and is not applicable to previous models which assume constant 
torque between stator and rotor. 


NFpR 


1 + 


NR-^Cs 


NFpR 


'}^W) 
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Conversely, at low loads (Cl/Cr < !)■ 


V'T/low — 


NFpR 


, jV-R'Cs \ 

^ (Ci-+C-r)^p J 


NFpR 


^1 , nr^Cs \ ■ 

. + Cnll J 


( 33 ) 


The torque measured at high speeds Is approximately 20 pN-nm |[T5]. Then, from Equation ||^, the non- 
dlmenslonal number NR'^Cs/Cr^I ~ 10. Also, Equation shows that the torque and speed are not linearly 
dependent on the number of stators at low loads. This Is consistent with previous experimental observations Il34l . 

The above approximations to the full model are primarily laid out for expository purposes, to Introduce general 
properties of the model. As such, there are several limitations due to the assumptions made, which we outline 
below. 


Firstly, the observations from Sowa et at. are under question given that motors with more than one (or a few) 
stators may not have been considered In low-load measurements. In order to determine whether the stators do 
Indeed act Independently or in synchrony, experiments that directly account for the number of active stators must 
be performed. The extension of our model to N stators is dependent on the assumption that all stators step in 
synchrony, and the above results will not hold if stators are independent stochastic steppers. Because of the 
recent measurements made on single-stator motors, we have chosen to largely focus on the explicit modeling of 
the intrinsic mechanism of torque-generation in the BFM, rather than on the potential interaction between stators. 
While a thorough exploration of the dynamics of multiple-stator motors is extremely worthwhile, it is out of the 
scope of the current paper. 

There are also some inconsistencies in the approximation of an infinitely stiff spring at intermediate and low loads 
when the motor is not in the mechanically-limited regime. We have that « 7 nm, i? « 20 nm, and the non- 
dimensional number NR^Cs/Cr^I ~ 10 at low loads. Then from Equation |^, the ratio of the stator and rotor 
drags Cs/Ci? ~ 1- Using this, we can estimate the viscosities of the stator and rotor. The drag of the rotor is given 
by 
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Cr = -^VrR^ , (34) 

where rju is the viscosity of the rotor. . An estimate for the drag coefficient of the stator can be obtained by using 
Equation and the maximum observed speed of the wild-type motor (200 Hz), yielding Cs = 0.017 pN-nm- 
s/rad. Since the drag ratio of the stator and rotor is approximately unity, Cr ~ 0.017 pN-nm-s/rad. Using j^, the 
viscosity of the rotor is approximately 2 Poise. This is consistent with the fact that the cytoplasm is a mixture of 
water and proteins (T]. 


Likewise, the drag coefficient of the stator loop with a lever arm of length £p is given by 


TT risip 

3 log(;^) - 0.66’ 


(35) 


where ijs is the viscosity of the stator and Rp is the radius of the stator loop. Then, using and a length to 
width ratio £p/Rp « 10, we estimate the viscosity of the stator to be 800 Poise. This value is three orders of 
magnitude higher than the viscosity of a regular lipid membrane and two orders of magnitude higher than the 
viscosity of a biological membrane. 


Moreover, the torque-speed curves produced by the stiff-spring approximation are linear, in contrast to the 
concave-down torque speed curves observed in single stator motors |T5]. Since the IMF enters only through 
the dependence of the torque, the speed-IMF curves are also linear at all loads, again in contrast to the recent 
100 nm bead experiments on a single stator fTSl . 
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1.1.3 Addition of chemicai kinetics 


The linear torque-speed curves of a mechanically rate-limiting model elucidate the importance of the inclusion of 
ion-binding kinetics. In this case, these are events related to the binding of a cation from the periplasm to Asp32 
and the unbinding of the cation from Asp32 into the cytoplasm. In this section, we recompute model torque- 
speed curves explicitly including the dwell times corresponding to the ion-binding and unbinding events between 
the power strokes. 


As done by Meacci and Tu [lO], a torque-generation cycle is divided into two parts: (i) moving time Tm and (ii) 
waiting time T^,. Assuming that the ions bind only when the stators are around the minimum of the respective 
free energy potentials, we may use the above model during Tm and sample from an exponential distribution 
at the end of each moving step. 


In a motor with a single stator {N = 1), the instantaneous torque is obtained from Equation as 


FpR 


(1 + 


R^Cs ' 

(Ci+CH)^p ' 

AG 


(36) 


The force applied by the proline hinge is given as Fp = = - 


27r 


Acj) 


2q IMF 
A(j) 


The time required to move an angular distance of — in step i can be calculated from Equation 
relationship w = r/(Ci? -f Cl)- 


where q is the charge of the ion. 

through the 


= T = — 

m m 26 


(Ci? + Cl)( 


1 + 


RHs \ 1 

{Ql F U)Pp) FpR- 


(37) 


Let TC, be the waiting time that follows a step i. During the waiting time, the instantaneous torque is zero. When 
the system reaches the steady state, the average torque (r) can be obtained via a time average 


(t) = lim — 
' ' T^oo T 


-(i')di', 


where T is large. If there are N steps in time T, then there are N waiting times. Therefore, Equation 
approximated as 


(38) 
can be 


(t) = lim — 
' ' T^oo T 


r(i')dt' 


Yi:=rrT„ 


T 


i—1 ~ A^i—1 


(39) 


tT„. 


T 

rr 


1_ rp 

N ^i—1 ^ I 


The average dwell time can be approximated as (T.^,) = ^fcoexp 

proportionality constant related to the rate of hopping of the ions. Using this. Equation 

(r) = 


tT„ 


, where kc is a 
reduces to 

(40) 


Tm + (Tju) 

Likewise, the speed of the rotor (or the bead in the large spring constant limit) can be written as a time average 


' dt ' 


lim — 
T-s-oo T 


dt' 


2Tr 


dt' = 


T 

rr 


_ T 

26Tm ^ 

I y-vA/' rp. 


2tt 


T„ 


26Tp, ™ 

Tm + {Tw) 


(41) 
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Figure 10: Contact forces. (A) Case 1: a force F is applied to object A which is in contact with object B. (B) 
Free-body diagrams for A and B arranged as in Case 1. (C) Case 2: N objects of type A are in contact with 
object B, each with a force F applied to it. (D) Free-body diagrams for each object of type A as well as B for 
Case 2. 


where the instantaneous speed is || x during the power stroke and zero during the dwell time. 

The torque during a power stroke r, the time taken by a single step and the dwell time all depend on the 
IMF. Thus, and ||^ point to the existence of nonlinearity in the torque-speed and speed-pmf curves. These 
curves are presented and explained in the Results and Predictions section of the main text. 


2 Overview of steric forces 


In this section, we describe our modeling of the steric forces between the stator and the rotor, as well as a general 
discussion on the nature of contact/steric forces in a low Reynolds number environment. The behavior of objects 
moving at low Reynolds number is counterintuitive. When the Reynolds number is small, viscous forces dominate 
over inertial forces and inertia can be ignored In the following, we illustrate some of these properties using 
a simple linear momentum balance. We can then extend this analysis to angular momentum balances, which are 
directly relevant to the BFM. 

Consider a force F pushing an object A which is in contact with a larger object, B, as shown in Figure [T^ We 
denote the drag coefficients on objects A and B as (a and (b, respectively. Let us consider the following two 
situations. 


Case 1 : if is applied to A, how much force is transferred to B when they are in contact, as shown in FigurejT^? 
The corresponding free body diagrams for A and B are shown in Figure [To^ ; in the steady state, the force 
balances can be written as: 


F - Fc- Cava = 0 
Fc — Cbvb = 0 . 


(42) 


In Equation va and vb are the velocities of objects A and B, respectively, and Fc is the contact force 
between objects A and B. When the two objects are in contact and are moving together, the following ‘contact 
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condition’ ensures that the velocity of the objects are equal: 


( 43 ) 


VA = VB- 


In this case, the force transferred by A to i? (the contact force, Fc), is obtained by solving j4^ : 

F= ^ . 

1 + Ca/Cs 


(44) 


When the drag on B is large, i.e., (a/(b < 1, almost all of the force F is transferred to the object B. Conversely, 
when Ca/Cb > then most of the force is consumed to drag object A with little force transferred to object B. 


Case 2: when there are N objects of type A in contact with B, each with a force F applied to them. This case is 
shown in Figure [Top. As in the first case, we are concerned with how much of F is transferred to B. Again, using 
free-body diagrams (shown in Figure pTop), we can write the equations of motion for the objects as 

F-F^-Cav\ = 0 

^ . (45) 

X! - Cbvb = 0 , 


where F^ is the contact force between i*'* object of type A and object B. As before, when the objects are in 
contact, the contact condition ensures that v\ = vb for all i. Therefore, the force transferred to B by TV objects 
of type A can be derived from as 


F, 




NF 


I + NCa/Cb' 


(46) 


When Ca/Cb < 1, Fc « NF. Therefore, the force transferred is multiplied by the number of objects pushing B. 
However, when Ca/Cb > 1, then the force transferred is Fc « F/(Ca/Cb) « 0. Finally, when Ca/Cb « 1, then 
the force transferred is Fc « F (i.e. the force transferred in this system is approximately the same as a single 
object A pushing B). 


The above properties of contact forces can be applied to the BFM by identifying object A as the stator and B as 
the rotor-load system (i.e., the rotor and the bead), and F as the internal force generated by the proline hinge 
pushing the stator from its straight to bent state. 


For large loads (e.g., large beads), almost all of the force generated by the stator is transferred to the rotor; 

i.e., the torque is close to stall. Moreover, as more stators are recruited, the force transferred increases linearly 
with the number of active stators, in accordance with the observed linear speed dependence on the number of 
stators near stall torque. By contrast, at zero load, if there exists a situation where the stator and rotor drags are 
comparable, then the above analysis suggests that the force transferred during a single step is equivalent to a 
single stator pushing the rotor. This suggests that, if the assumption that stators step in synchrony holds, torque 
and speed at low loads may be independent of the number of active stators. 


In addition to the above, contact forces also have the following important implications for the BFM: 


1. Because it operates at low Reynolds number, the rotor moves only as long as it is pushed by the stator. 
This assures that the rotor never moves faster than the stator. 

2. Experiments based on torque-speed curves alone may never be able to detect the number of operating 
stators at different loads (e.g., if both torque and speed are independent of the number of stators at low 
load). Rather, one likely needs chemical markers such as GFP tags as used in (30] to identify the number 
of docked and engaged stators for a given load. 

3. The torque generated by the BFM depends on the bead size and is not constant across applied loads, as 
it was considered to be by previous models (331120] [191 . 
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Figure 11: A schematic of the stator potentials. 


3 Numerical implementation 


In the above sections, we have provided the equations of our model, as well as analytic solutions for reduced 
deterministic approximations. In this section, we provide details on numerical implementation of the full stochastic 
model, utilizing several arguments from our discussion on steric forces. 


3.1 Interaction potentials 

In this section, we provide explicit forms of the free-energy potentials in terms of the order parameters we choose 
to describe the motion of the stator and the rotor. 


3.1.1 Stator potentials 

The order parameter describing the motion of the stator is the angle subtended by the stator loops (j)^^ with 
respect to the vertical MotB ion channels. When two ions bind to the Asp32 residues, the two loops undergo a 
conformational change from their straight = 0°) to bent state = 20°). In this work, we have assumed 
that the stator loops move in-phase. Therefore, as previously, we model the two stator loop configurations using 
a single collective parameter (j)s. 

Before the ions bind to the Asp32 residues, the motion of the stators is governed by the potential Gi{4>s) where 
the minimum is around 0s = 0° as shown in Figure [il] When two ions bind to two Asp32s, the stator potential 
switches from Gi to G 2 . This compels the stator angle to move from 0s = 0° to 0s = 20°. 
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Figure 12: A Schematic of the Geometry of the Rotor-Stator Contact during the half step. The distance between 
the nearest FliG and the stator is then given by xrs = ^rs + RSr - £p4>s- The contact force is zero, whenever 
xrs > Xrs- Moreover, the stator never is ahead of the rotor. This is ensured by choosing a time step such that 
the fluctuations from the random noise terms are small. 


During this transformation, the loop pushes the rotor via a steric force. At the end of the conformational change, 
when the loops are at the minimum of the potential G 2 , the two ions bound to two Asp32s exit into the cytoplasm. 
The potential then switches back from to Gi and the loops traverse back to (j)s = 0°. During this time, the loops 
apply a contact force on the same FliG as in the previous sub step. 


As noted previously, thermal fluctuations are of the same order of magnitude as the free-energies considered 
and the precise form of the potentials is not important. In our simulations, we choose: 


and 


GM) 



if (/)< 0 

Fp(j) 

if 0 ^ (pruax 

FpCffrnax “F 4^max) 

if ^ > (jy^max 


G 2 W 


\ -Fpi^P) 

Fp(l)max “F 4^max') 


if (/) < 0 

if 0 ^ ^ (j^niax 

if 0 ^ 4^raax- 


3.1.2 Rotor-stator interaction potentiai 


The steric force between the stator and the rotor can be simulated 
distance Xrs. 


Vrs{<Psi (^r) 



{ROp—^rcps)^ 

Xrs 


using a soft linear repulsive force with a cutoff 


if 0 < a: < Xrs 

otherwise, 
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The torque on the rotor can be obtained as rcontact 


_ _ dVns ■ 

~ dBu ■ 


'r’contact(2^) — 


—FrsR 

0 


R9ji—lp4>s _ 

Xrs 


FrsR 


1 _ XRS A 

^ Xrs I 


if 0 < Xrs < Xrs 

otherwise. 


Likewise, the reaction torque on the stator is Frs is the maximum force that can be applied by the proline 

hinge and xrs = XRs+ROR-ipcjis denotes the distance between the position of the stator loop and the nearest 
FliG (see FigureQ^. All other parameters are defined as described previously. 

Note that we assume that weak electrostatic forces place the stator at most 0.5 nm from the nearest FliG on the 
rotor {Xrs ^ 0.5 nm) prior to the start of the power stroke, as shown in Figurejl^ During the half step, the stator 
moves from (/)s = 0° to (ps = 20°. In this process, a contact force is applied on the rotor when 0 < xrs < Xrs- 
The contact force is zero when xrs > Xrs- 

All previous theoretical studies of the BFM have chosen to model the interactions between the stator and the 
rotor as a load-independent force. In the above sections, we have shown that the dependence of steric forces on 
the load reproduces many of the mechanical characteristic features of the motor found in experiments (e.g., the 
nonlinearities in the relationship between speed and stator number at high loads). 

While the exact form of the potential used to model the rotor-stator interaction is not very important, we note 
linear potentials cannot be used. Such potentials result in a constant torque, independent of Cl, and thus does 
not reproduce several properties of steric forces. 


3.2 Kinetics of ion-binding 

In this section, we will describe how to model ion-binding events as part of the Langevin equation framework. 
In many problems related to motors, the ion-binding events may occur only in a certain window of a continuous 
coordinate describing the mechanical motion of the system. In the BFM, there are two main ion-binding events 
involved in the torque-generation cycle. 


1. Two cations (here, protons) from the periplasm bind to the Asp32 residue on two MotBs or unbind into the 
cytoplasm when the cytoplasmic loops are straight. 

Asp32- -t- ^ Asp32-H (47) 

2. Once the protons are bound, the cytoplasmic loops undergo conformational change from the straight state 
into the bent state. At the end of the conformational process, the two channels close with respect to the 
periplasm and instead open towards the cytoplasm. Once this occurs, the two protons unbind from Asp32 
residues into the cytoplasm (or, in the reverse reaction, bind from cytoplasm to Asp32). 

Asp32-H 4^ Asp32- -f H+(48) 

ki2 

First, we will describe the case for the kinetics at equilibrium (i.e., at zero IMF), followed by the procedure to 
model the kinetics under a non-zero IMF. 
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3.2.1 Equilibrium kinetics under zero IMF 


In the following we consider the proton-driven motor of E. coli, and so the IMF concerned is the proton motive 
force (PMF). When there is no PMF, the forward and backward reaction rates for the reaction in j47) must satisfy 


ki2 

k2i 


Let the equilibrium dissociation constant of reaction 


[Asp32-H] 
[Asp32] ■ 

be defined as 


K^ = 


[Asp32][H+ 


periplasmJ 


[Asp32-H] 


(49) 


(50) 


If pK^ = - logiQ K^, then the forward and backward rates should satisfy 


— = lo(P'^a-PHperiplasm) ^ 

fel 


(51) 


Note that the acid dissociation constant value of pK^ should be determined from experiments. A similar relation 
to the above can be derived for the reaction in 


3.2.2 Kinetics under non-zero IMF 

To satisfy detailed balance when the PMF is non-zero, the kinetic coefficients for the reaction should satisfy 

hi = 10(PKr-pHpePp,a3J exp ( ^^P + Gi-G2 \ ^g^) 

^21 V ksT J 

There exist multiple choices for the expressions ki 2 and k 2 i. For convenience and symmetry, we choose the 
following for the rate constants : 

ki2 = exp ^ 

k2i = exp (^{l + a)(^ 


4 Electrostatic steering 


In this section, we provide more detailed calculations for our electrostatic steering hypothesis, as shown in the 
main text. We also provide support for the dipole approximation via comparison with an analogous calculation 
with point charges. 

We note that, because no structure of the stator is yet available, the purpose of these calculations is quite 
qualitative. Our goal is to simply to predict what the electrostatic energy landscape should look like in order to 
support our steering hypothesis. 
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4.1 Dipole approximation 


As a first-order approximation, we consider the relevant charges on the FliG helix and the stator loops to be 
dipoles. We denote the rotor dipole moment as where fc G {1... 26} enumerates the number of FliGs on 
the rotor periphery. Likewise, the stator dipole moment is denoted as Pg^ , where N enumerates the number 
of stator units, and i the number of loops on a single stator. Here for ease of exposition, we show electrostatic 
calculations for a single stator with a single loop (i = 1 and N = 1), but note that this calculation can be easily 
extended to the full stator model. 

We calculate the electric field felt by the stator loop as: 

26 

E = Y,Ek, (55) 

k^l 

where 

^ 4^e|rfc|3 

Here e is relative permittivity of cytoplasm and fk is a vector quantity which denotes the distance between the 
stator loop and the fcth FliG. 

From Ek, we calculate the interaction energy between the dipole on the stator loop pg and the fcth FliG pj^ as; 

Uk = -ps ■ Ek (57) 

= [(Pr ■ Ps) - Hpr ■ ^kWs ■ fk)] . 

Similar to the calculation for total energy, we have 

26 

t^tator = 'y ' Uk- (58) 

k=l 

Note that, as the distance from the stator loop increases, the terms in the total energy sum drop off as and 
so the contribution by FliGs located far from the stator loop is not appreciable. 

To calculate the total energy in the system, we also add in the interaction energies between pairs of FliG 
molecules. However, we note that their relative positions do not change as the stator rotates, and so this consid¬ 
eration results simply in a translation of the entire landscape and has no effect on the topology. 

[/total = [/stator + H A -^-13 \^PR ' PR^ ~ ^^PR ' E,j){Pr ' fi,j)] , (59) 

where is the displacement vector between the ith and jth FliG. 


4.2 Comparison with calculation using point charges 

We also performed the above calculations using point charges for the relevant residues on FliG and MotA loops. 
Because of the uncertainty in the position of the charges, we are interested primarily in a qualitative affirmation 
of our approximation — that is, the existence of a gently sloping, relatively wide energy well. 

As in our dipole calculations, we consider the charged residues on the FliG to be positioned at | to the horizontal 
(i.e., the charges are positioned along the dipole as shown in Figure 3 in the Main Text). Similarly, we position 


34 




0R (rad) (rad) 



Og (rad) 


Figure 13: Energy landscape during ‘electrostatic steering’ calculated using point charges rather than dipoles. 
All parameters are defined as in the main text; the stator is centered at (21,-2,1) and the rotor is centered at the 
origin (all distances in nm). 



the charges on the stator loop along the dipole, positioned radially outward from the stator center. Given this 
configuration, we calculate the electrostatic energy as 
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( 60 ) 


We screen charges using a Debye length of Ad = 0.5 nm. As before, denotes the displacement vector 
between charge i and charge j. Note that, similar to our dipole calculations, the first two terms of the energy are 
invariant to rotation of the stator. Therefore, they only serve to translate the entire energy landscape and do not 
affect the topology. 

As shown in Figure the topology generated is indeed similar to that using the dipole approximation in the 
main text. The removal of certain charges may modify the energy landscape via a flattening or widening of the 
energy well. In conjugation our electrostatic steering hypothesis, this may point to a possible explanation for 
experimental studies which show that mutation of charged residues reduce, but do not eliminate, motor function. 
Our methods can easily be extended using positional information from solved structures to more quantitatively 
analyze the effects of these charges. While such an investigation is certainly warranted, it is out of the scope of 
the current paper. 


4.3 Centering the stator 

There is little structural information for the stator of the BFM. Therefore, in our calculations, we choose a possible 
position of the stator which is consistent with our model and experimental results. Firstly, we center the stator so 
that all four stator loops are able to access the FUG molecules. The occasionally observed backsteps imply that 
both pairs of loops must be able to execute a power stroke during motor function. 

Preliminary calculations strongly suggest that the loop charges cannot be in plane with the charges on the FUG. 
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The close proximity of the charges results in an electrostatic interaction far too strong to be feasible given the 
known efficiency of the BFM (i.e., far too much energy would be wasted in “letting go”). In-plane configurations 
which do not present both extremely steep and deep wells must place the stator out of reach of the FliGs. 

Given these considerations, and noting that the stator and rotor have radii of respective lengths --^20 nm and '■•^2 
nm, we center the stator at {xs, ys, zs) = ( 21 , - 2 , 1 ) using a coordinate system with the rotor centered at the 
origin. 


5 An Ising model to describe conformational changes in FliG 


The existence of motor backstops in the absence of CheY-P has been largely attributed to “microscopic reversibil¬ 
ity”. Given this constraint, there are three possibilities in explaining the backstop; (i) A backward transit of an ion 
from the cytoplasm to the periplasm, (ii) A sufficiently strained hook relaxes by moving the rotor backwards, (iii) 
A fluctuation which changes the orientation of FliG from its CCW to CW orientation. 

The first two possibilities contradict our proposed model. For example, if the first case were to be correct, then 
a backstep would correspond to a reversal of the entire conformational change process. In this case, an ion 
bound to the Asp32 from the cytoplasm should restore the MotA helix from its bent state to the straight state in 
the presence of the ion. This is in direct contradiction with one of our model assumptions that MotA relaxes to 
the bent state due to the rearrangement in the hydrogen bonds caused by the binding of the ion to the Asp32 
residue. 

The second possibility (that the relaxation of an elastically strained hook due to several sequential forward steps 
leads to an occasional backstep) is also not feasible within the context of our model; a contact force guarantees 
that the rotor always follows the stator. 

In this work, we attribute the molecular basis for the backstep to the third possible scenario; a conformational 
change in one (or more) individual FliGs on the periphery of the rotor. A FliG can exist in two states; in our model, 
these states correspond to two orientations of the FliG dipole vector. Our electro-steric model, in conjunction 
with an Ising model corresponding to the states of FliG, explains the existence of occasional backstops in a 
self-consistent manner. 

Briefly, a backstep results from a fluctuation of FliG from the CCW to the CW position. In this case, when the 
motor is predominantly moving in the CCW direction and whenever a FliG changes its state from CCW state 
to CW state and is close to a stator, then the stator—using MotA loops 2 and 4—applies the contact force and 
pushes the FliG in the backward direction. 

While the configuration of the FliGs in the CW orientation is fairly well agreed upon, the exact CCW configuration 
is still under debate. At least three possible orientations have been suggested, including directions orthogonal 
or 180° with respect to the CW orientation. Despite this uncertainty, it is generally believed that there are two 
significantly different orientations for each FliG. 

As in other models, we describe the transitions between the two FliG orientations by a one-dimensional periodic 
Ising model consisting of 26 spins (corresponding to 26 FliGs). Each individual spin can exist in two possible 
states corresponding to the two orientations of FliGs; Si = -1-1 (CCW) and Si = -1 (CW). The Hamiltonian is 

N N 

H = —J SiSj — h Si (61) 

i,j — l i—1 

where J denotes the nearest neighbor pairwise interaction energy and h denotes the field biasing the FliGs to 
preferentially orient in a certain direction. Let {s^} = (si, S 2 , S 3 ,..., S 26 ) denote a possible state of the rotor. The 
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probability of such a state is given by the Boltzmann distribution 




where Z{J,H,P) is the partition 


function. The partition function for a one-dimensional Ising model is obtained exactly using the transfer- 
matrix approach, and is given by 

Z{J,h,p) = X^ + X^ 


where 


X± = cosh(/3h) ± sinh^(/3/i) -|- . 


(62) 

(63) 


In order to understand the back steps while the motor is running predominantly in CCW direction, we calculate 
the probability that any FliG (or any spin) points in the CW direction. This is given by 

p-/ 3 Jr((si,...,Si_i,-l,Si + i,...,S26)) 

p(s^ = -1) = - (64) 

where the summation is over all possible states fixing Si = —1. In the following, we derive an analytical expression 
for the probability p(si = -1). 


Let us denote the probability for the spin s* = +1 as p(si = -l-l). Given the probabilities p(si = -l-l) and 
p(si = -1), the average value (si) of the spin Si can be calculated as 

(si) = (l)p(si = -hi) -L (-l)p(si = -1) = 1 - 2p(si = -1). (65) 

Moreover, the derivative of the partition function Z with respect to the field h yields 

Nf3Z dh Z ^ ^ 

{Si} * 

Therefore, using and |^, the probability that spin = -1 can be obtained as 

1 _ 1 dZ 

, X ^ JV/3Z dh 

p{si = -1) = - ^ - . (67) 

The derivative of the partition function with respect to the field h can be evaluated using the derivatives of the A± 
functions, which are given by 


9A± 


/?sinh(/3/i) [ ± 


cosh(/3/i) 


\ sinh^ [fih) + ) 


( 68 ) 


Within the context of our model, the likelihood of a flip in FliG conformation corresponds directly to the likelihood of 
a backstop. That is, the probability that a FliG is in the CW state is given by |[^. This is because the fundamental 
mechanics of a backstop is the same as that of a forward step; the difference between these two scenarios is 
wholly described by the difference in FliG configuration. 

Suppose that the motor is moving primarily in the CCW direction. This means that the FliGs in close proximity 
to the stator loops are oriented to favor CCW rotation (that is, to favor interaction with MotA loops 1 and 3). 
However, if a FliG close to the stator is oriented to favor CW rotation, then the FliG interacts with MotA loops 
2 and 4, resulting in a step in the CW direction (a “backstop” when the motor is moving primarily CCW). This 
probability is given by (Ft) . 

For the choice of the energy scale /3 J = 2 at room temperature and a biasing field of ph = 0.05, the probability 
that any FliG is in the CW state is p{si = -1) = 0.08. That is, on average, 8 out of every 100 torque-generating 
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Figure 14: Probabilities of backward steps relative to forward steps can be estimated using statistics from stepping 
experiments (data shown from |^). The area under the curve for steps below and above 0° corresponds to 
forward and backward step probabilities, respectively. The areas on either side are calculated as highlighted 
in gray. Curves are truncated where each shows a second “peak”, corresponding to two steps being blurred 
together due to limitations in experimental resolution. The relative areas calculated for the data above suggest 
that the probability of a backwards step is 7.3%. 


cycles will result in a backstep. Areas under the curve corresponding to forward and backward steps from data 
collected by Sowa etal. |29] indicate thatp(si = -1) ~ 0.073 (Figure [14 '. 


The above analysis provides an explanation for how backstops in the absence of CheY-P can arise from fluctu¬ 
ations in FUG configurations. However, it does not take into account how the above probabilities are effected by 
load or PMF. For example, the timescale of a single step depends on the load. If one assumes that this step is a 
backstep, then this particular FUG is pinned in the “backwards” orientation for the duration of that step. This CW 
defect in the FUG ring can affect the switching probabilities of the neighboring FliGs, resulting in further defects 
along the ring. Therefore, at higher loads, the probability of two (or more) subsequent backsteps may not be 
negligible. However, a complete understanding of the above requires a far more detailed analysis of an Ising 
model in conjunction with the proposed electrosteric model than is within the scope of this work. 
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